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X-autosome translocation (XY Y,) has been reported in distinct groups of vertebrates suggesting 

that the rise of a multiple sex system within a species may act as a reproductive barrier and lead to 
speciation. The viability of this system has been linked with repetitive sequences located between 

sex and autosomal portions of the translocation. Herein, we investigate Oecomys auyantepui, 

using chromosome banding and Fluorescence In Situ Hybridization with telomeric and Hylaeamys 
megacephalus whole-chromosome probes, and phylogenetic reconstruction using mtDNA and nuDNA 
sequences. We describe an amended karyotype for O. auyantepui (2n = 649 654/FNa = 84) and report 
for the first time a multiple sex system (XX/XY,Y,) in Oryzomyini rodents. Molecular data recovered 

O. auyantepui as a monophyletic taxon with high support and cytogenetic data indicate that O. 
auyantepui may exist in two lineages recognized by distinct sex systems. The Neo-X exhibits repetitive 
sequences located between sex and autosomal portions, which would act as a boundary between 
these two segments. The G-banding comparisons of the Neo-X chromosomes of other Sigmodontinae 
taxa revealed a similar banding pattern, suggesting that the autosomal segment in the Neo-X can be 
shared among the Sigmodontinae lineages with a XY,Y, sex system. 


Abbreviations 
mtDNA Mitochondrial DNA 
nuDNA Nuclear DNA 


Cytb Cytochrome b 

Col Cytochrome C Oxidase Subunit I 
FGB-I7 _ Beta-fibrinogen intron 7 

BI Bayesian Inference 

ML Maximum Likelihood 

FISH Fluorescence In Situ Hybridization 
ITS Interstitial telomeric sequence 

2n Diploid number 

FNa Autosomal fundamental number 
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CH Constitutive heterochromatin 
het-ITS Heterochromatic-ITS 
APRT Adenosine phosphoribosyltransferase 


HSA Human whole chromosome probes 
MMU Mus musculus 

HME Hylaeamys megacephalus 

OAU Oecomys auyantepui 


OCA-PA Oecomys catherinae From Para 

OCA-RJ — Oecomys catherinae From Rio de Janeiro 
OPA-A Oecomys paricola Cytotype A 

OPA-B Oecomys paricola Cytotype B 

OPA-C Oecomys paricola Cytotype C 


CLA Cerradomys langguthi 
NVO Neacomys vossi 

NEL Neacomys elieceri 
NXI Neacomys xingu 
NMA Neacomys marajoara 
NPA Neacomys paracou 
NSP-E Neacomys Sp. E 

NAM Neacomys amoenus 
TNI Thaptomys nigrita 
AMO Akodon montensis 
ASP Akodon Sp. 

NLA Necromys lasiurus 
OAM Oxymycterus amazonicus 
BBR Blarinomys breviceps 


Chromosomal rearrangements are drivers in karyotypic evolution and are often associated with speciation’. 
Mammals are known to exhibit a stable sex determination system, but distinct sex-autosome translocations 
may have triggered the separation of Theria and Prototheria (monotremes) (190 MYA) and between Euthe- 
ria (placental mammals) and Metatheria (marsupials) (166 MYA)°. Although the euchromatic region of the 
X chromosome is considered conserved among highly rearranged karyotypes of placental mammals’, recent 
investigations on Arvicolinae (Myomorpha) rodents have shown that the X chromosome has undergone several 
intrachromosomal rearrangements, such as centromere shifts, peri- and paracentric inversions, that were also 
accompanied by repetitive sequences’. Regardless of whether chromosomal rearrangements are the primary cause 
of speciation”*, or whether karyotypic divergence between closely related species are a casualty of the speciation 
process”’°, the most deleterious among the speciation-linked rearrangements’? are tandem translocations, 
reciprocal translocations'*"* and X-autosome translocations’*"*. 

The rise of an X-autosome translocation is subordinated to the same epigenetic mechanism that guarantees 
dosage compensation between normal females (XX) and males (XY) by silencing one of the Xs in females’’. In 
this type of event, the inactivation progress in one of the X chromosomes of females!® spreads to the autosomal 
segment translocated to the X, silencing genes in the autosomal portion’”’ generating deletion/duplications with 
deleterious effects”. 

Although deleterious effects of sex-autosome translocations have been described in the literature for humans 
and mice (e.g., male sterility; embryonic lethality)'**!”’, this type of chromosomal rearrangement has been 
reported in natural populations of distinct groups of vertebrates, such as fish?*”*, anurans* and mammals’””*. 
The presence of intercalary heterochromatic blocks between autosomal and ancestral X chromosome segments 
could suppress the X-inactivation progress in the autosomal segment, allowing viability in this system!°?”-. 
Several studies have shown the presence of heterochromatic blocks, telomeric repeats and/or rDNA (ribosomal 
DNA) clusters in different mammalian lineages that exhibit X-autosome translocation, for example in bats (Chi- 
roptera) of genera Artibeus, Carollia, and Uroderma”-*; in rodents (Rodentia) of genera Nannomys (Muridae)'”, 
Proechimys (Echimyidae)**~*°, and Taterillus (Muridae)'*; in diprotodont marsupials (Diprotodontia) of genus 
Wallabia*®; and in ruminants (Artiodactyla) of genera Antilope and Gazella’’. 

In rodents from the Brazilian Amazon, the XX/XY,Y, multiple sex system has been reported only in two 
genera from the Echimyidae family: Lonchothrix** and Proechimys**-*. In Lonchothrix emiliae, the multiple sex 
system was identified based on classic banding**, while in the Proechimys taxa it was detected by FISH (Fluo- 
rescence In Situ Hybridization) with whole chromosome probes (chromosome painting) from P. roberti and P. 
goeldii**. In Sigmodontinae rodents (Rodentia, Cricetidae), the Oryzomyini tribe currently comprises 29 genera 
and is the most diverse of the 11 tribes within the subfamily*’', but multiple sex systems are acknowledged 
solely in representatives of the Akodontini, Phyllotini and Reithrodontini tribes: Deltamys kempi (Akodontini) 
exhibits a X,X,X,X,/X,X,Y sex system due to a translocation involving chromosomes 2 and Y”; Salinomys deli- 
catus (Phyllotini) shows a XY,Y, system“; and Reithrodon (Reithrodontini) exhibits a XY, Y, system (Uruguay 
population) and a Neo-XY system (Brazil population)“. 

In Oryzomyini, the genus Oecomys has been particularly challenging in taxonomy, distribution patterns and 
speciation mechanisms. Comprising 19 species to date, Oecomys has been investigated using several approaches, 
such as morphology, nuclear DNA (nuDNA), mitochondrial DNA (mtDNA), and cytogenetics, which have 
shown that some lineages correspond to species complexes*”**>!. Oecomys auyantepui has been recognized as 
a monophyletic lineage and a monotypic taxon***’. The species is distributed from southeastern Venezuela to 
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Oecomys auyantepui (OAU, 2n=653/FNa=84) 
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Figure 1. Oecomys auyantepui (2n=653/FNa = 84) (a) G-banded karyotype with chromosome painting 
revealed by Hylaeamys megacephalus (HME) whole chromosome probes™, and (b) C-banded karyotype. An 
asterisk indicates a centromere. 


north-central Brazil, in the Guiana subregion of Amazonia”’, and exhibits two sympatric populations with distinct 
diploid numbers (2n) of 64 and 66 and autosomal fundamental numbers (FNa) of 110 and 114, respectively’. 
A third karyotype of 2n =72/FNa = 80 was described**. In addition, an interstitial telomeric sequence (ITS) was 
identified at the centromeric region of the bi-armed X chromosomes in karyotypes with 2n = 64 and 66, which 
suggests that chromosomal rearrangements have driven the evolution of this chromosome in O. auyantepui”. 

It is noteworthy that cytogenetics studies with Oecomys have shown a substantial diversity in 2n and FNa, 
ranging from 54 to 86 and from 62 to 140, respectively*?178"°, However, specific events that shaped extant 
karyotypes remain unclear for most species, except for O. catherinae from Para (OCA-PA; 2n=62/FNa =62), O. 
catherinae from Rio de Janeiro (OCA-RJ; 2n = 60/FNa = 62), O. paricola cytotype A (OPA-A; 2n =72/FNa=75), 
O. paricola cytotype B (OPA-B; 2n=70/FNa=75), and O. paricola cytotype C (OPA-C; 2n=70/FNa=72) that 
were investigated by chromosome painting with Hylaeamys megacephalus whole chromosome probes (HME; 
Oryzomyini)*”*’. In addition to elucidating the chromosomal rearrangements that occurred in these species, 
the chromosome painting analysis helped to delineate taxonomic limits, as the authors*”*° were able to identify 
a hidden diversity and proposed that O. catherinae and O. paricola “eastern clade” were composed of two and 
three species, respectively. 

Considering the evolutionary force of chromosomal rearrangements regarding speciation and diversification 
of species, we set out to investigate if the emergence of a new sex-system triggered the speciation process in the 
monotypic taxon Oecomys auyantepui. 

In order to achieve this goal, we used classic cytogenetics, telomeric and HME whole chromosome probes”, 
mtDNA (mitochondrial DNA) and nuDNA (nuclear DNA) sequences. Here we discuss the chromosomal evo- 
lution of the genus, and report for the first time a multiple sex system (XX/XY,Y,) in Oryzomyini rodents. We 
also compared the taxa from the present study with other species analyzed elsewhere using the same set of 
probes!”°54-59, 


Results 

Classic and molecular cytogenetics. Oecomys auyantepui (OAU) has a 2n=649654/FNa = 84 karyo- 
type, with a multiple sex system (XX/XY,Y,). The autosomal set consists of 20 acrocentric pairs (1-20) and 11 
meta/submetacentric pairs (21-31). In females sex chromosomes were recognized as two medium-sized sub- 
metacentric Neo-X chromosomes; in males sex chromosomes were identified as one Neo-X and two Ys: Y, chro- 
mosome was a medium submetacentric (original Y) and Y, was a small acrocentric (Xp homologue) (Fig. 1a). 
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21 8q dist 
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Table 1. FISH results for Oecomys auyantepui (OAU; 2n=654/FNa = 84), as assessed based on hybridization 
with Hylaeamys megacephalus (HME) whole-chromosome probes™. 


The constitutive heterochromatin (CH) is distributed in the centromeric regions of almost all autosomes, the 
Neo-X and Y, chromosomes. The CH is a small region in most of the autosomes and the Y, chromosome has a 
large heterochromatic block in the long arm (Fig. 1b). 

Cross-species FISH with HME probes yielded 42 signals on the OAU chromosomes (Fig. 1a, Table 1, see 
Supplementary Figs. 1 and 2). Ten autosomal probes are conserved; of them, four (HME 8, 15, 24 and 25) 
hybridize to whole chromosomes of OAU (1, 10, 11 and 28, respectively) and six (HME 7, 12, 18, 20, 21 and 
26) hybridized with portions of chromosomes of OAU (21q, 23p, 25q, 25p proximal, 8q distal and 25p distal, 
respectively). Twelve autosomal probes show multiple signals in OAU: HME 1 hybridize to OAU 2, 19 and 23q; 
HME 2 hybridize to OAU 4 and 22q; HME 4 hybridize to OAU 14 and 24; HME 5 hybridize to OAU 7, 18 and 
31; HME 6 hybridize to OAU 21p and 30; HME (9,10) hybridize to OAU 5, 22p and 26p; HME 11 hybridize to 
OAU 9 and 20; HME (13,22) hybridize to OAU 8q proximal, 17 and 26q; HME 14 hybridize to OAU 13q distal 
and 15; HME (16,17) hybridize to OAU 6 and 27; HME 19 hybridize to OAU 12 and 13q proximal; HME 23 
hybridize to OAU 16 and 29. 

The HME 3 probe hybridizes to the short arm of the Neo-X chromosome (OAU Xp), and also hybridizes to 
OAU 3, and Y,; the HME X chromosome hybridizes to the long arm of the Neo-X (OAU Xq). 

Seven OAU autosomal pairs show hybridization signals to multiple HME probes: OAU 8 (HME (13,22)/21); 
OAU 13 (HME 19/14); OAU 21 (HME 6/7); OAU 22 (HME (9,10)/2); OAU 23 (HME 12/1); OAU 25 (HME 
26/20/18); OAU 26 (HME (9,10)/(13,22)) (Fig. 2). 

FISH with telomeric probes showed hybridization signals at the distal regions of all chromosomes, plus a large 
interstitial telomeric sequence (ITS) at the centromere of the Neo-X chromosome (Fig. 3). 


Phylogenetic analysis. A more detailed phylogenetic analysis of the genus Oecomys was already 
proposed“. Thus, in this work we focused on O. auyantepui and representatives of each Oecomys species/clade 
recognized in the literature***”>! (Supplementary Table 1). The genus Oecomys was recovered as monophyletic 
in the topologies obtained with the Cytochrome b (Cytb) dataset, the Cytochrome C Oxidase Subunit I (COT) 
dataset, and with the concatenated dataset (Cytb + beta-fibrinogen intron 7 [FGB-I7]), with high support values 
recorded only in the Bayesian Inference (BI) analyses (Figs. 4, 5, 6). In the Cytb topology, lineages of O. bicolor 
and O. cleberi were not recovered as reciprocally monophyletic, as well as lineages of O. mamorae and O. fran- 
ciscorum (Fig. 5). 

We found mean interspecific p-distances ranging from 6.47% (between O. concolor and Oecomys sp.2) to 
14.80% (between O. franciscorum/O. mamorae and O. matogrossensis), and mean intraspecific p-distances vary- 
ing from zero in several species to 6.2% in O. bicolor/O. cleberi (2.15% in O. auyantepui, particularly; Table 2). 
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Figure 2. FISH results obtained from O. auyantepui (OAU; 2n= 649654/FNa = 84), using HME whole 
chromosome probes”. Each box corresponds to an OAU pair shown in Fig. 1a that corresponded to more than 
one HME homologue. Single or multiple images are presented to exhibit full coverage with HME probes on 
OAU chromosomes. OAU chromosomal pairs identification are shown below the chromosomes, while HME 
probes are shown beside the chromosomes. An asterisk indicates a centromere. HME whole chromosome 
probes are shown in green (FITC) and red (CY3); the counterstaining is blue (DAPI). 





Figure 3. FISH results obtained from O. auyantepui (OAU; 2n=654/FNa= 84), using telomeric probes. Arrow 
indicates the interstitial telomeric sequence (ITS) at the Neo-X chromosome. Telomeric probes are shown in 
green (FITC); the counterstaining is blue (DAPI). 


In all three datasets, O. auyantepui was recovered as a monophyletic taxon with high support. The COI 
phylogeny was the only one that included all O. auyantepui karyotyped samples from this work and from the 
literature°'”’, as well as most sequences available on GenBank (Supplementary Table 1). In the COI topology, 
specimens of O. auyantepui formed a polytomy, with no resolution among most specimens, including those 
with similar karyotypes (Fig. 4). In the Cytb topology, the specimen N228 was recovered as the most divergent 
within O. auyantepui, with 5.20% of mean genetic divergence from its conspecifics. The remaining specimens 
formed a subclade with neither resolution nor support among most of the specimens. As in the COI topology, 
specimens with similar karyotypes did not nest in a subclade (Fig. 5). Finally, in the concatenated data topology, 
the specimen ROM 114,316 was the first one to diverge, and the specimen ROM 114,059 was recovered as sister 
to the 2n=654/FNa = 84 specimens included in this analysis. Although these latter specimens appeared as a 
subclade, there was no support for that (Fig. 6). 
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Figure 4. Bayesian Inference topology based on Cytochrome C Oxidase Subunit I. The numbers above 
branches indicate posterior probability values for Bayesian Inference analysis (only values > 0.80 are shown) and 
bootstrap values for Maximum Likelihood analysis (only values > 65% are shown). Bold numbers indicate the 
samples from this study. Sample data are provided in Supplementary Table 1. 


Discussion 
Chromosomal evolution and signatures in Oecomys (Rodentia, Sigmodontinae). As men- 
tioned above, there is a large variation in 2n from 54 to 86 and in FNa from 62 to 140 among the Oecomys species, 
with karyotypes mainly composed of one-armed chromosomes*”*”48592°306!, The variation in 2n and FNa 
occur both within and between species, indicating that fusions/fissions, pericentric inversions (or centromeric 
repositioning), translocations, and addition/deletion of constitutive heterochromatin are the main forces acting 
in the chromosomal evolution of this group of Sigmodontinae rodents. Thus, we used chromosome painting to 
make a comparative analysis of taxa in the present study (O. auyantepui) and in the literature (O. catherinae and 
O. paricola) to precisely identify the rearrangements among them*””°. 

As a result, chromosome painting analysis in the karyotypes of O. auyantepui (OAU), O. paricola (OPA- 
A, OPA-B and OPA-C) and O. catherinae (OCA-PA and OCA-RJ) (Supplementary Table 2) showed that the 
chromosomal variation in 2n from 60 to 72 and in FNa from 62 to 84 are due to 23 fusion/fission events, four 
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Figure 5. Bayesian Inference topology based on Cytochrome b. The numbers above branches indicate posterior 
probability values for Bayesian Inference analysis (only values > 0.60 are shown) and bootstrap values for 
Maximum Likelihood analysis (only values > 60% are shown). Bold numbers indicate the samples from this 
study. Percentage values are mean genetic distances (p-distances) between high supported selected lineages of O. 
auyantepui. Representatives of different lineages within O. bicolor, O. cleberi, O. roberti, O. mamorae, O. paricola, 
and O. catherinae complexes are indicated by the letters C (central clade), E (eastern clade), N (northern clade), 
NW (northwestern clade), S (southern clade), W (western clade), and Wm (westernmost clade), following’. 
Subscribed letters b and c indicate lineages attributed by** to O. bicolor and O. cleberi, respectively. Subscribed 
letters f and m indicate lineages attributed by** to O. franciscorum and O. mamorae, respectively. Cross (*) 
denotes karyotype information not obtained from the specimen included in the phylogenetic analysis. Sample 
data are provided in Supplementary Table 1. 


translocations, seven pericentric inversions and amplification/deletion of constitutive heterochromatin on two 
autosomal syntenic blocks plus the X chromosome (Supplementary Fig. 3), with only seven syntenic blocks con- 
served without detectable rearrangements. Remarkably, we observed that the rearrangements that differentiate 
OPA cytotypes (OPA-A, OPA-B, and OPA-C) from each other are different from those responsible for the vari- 
ability between OCA cytotypes (OCA-PA and OCA-RJ) and OAU. This suggests that the rearrangements mainly 
occurred in distinct syntenic blocks among these species (Supplementary Fig. 3). Consequently, we propose that 
each of these three species has evolved independently and has not followed the same path of rearrangements or 
the same chromosomes. 

Moreover, by detecting an elevated number of chromosomal rearrangements among three taxa (O. auyantepui, 
O. catherinae, and O. paricola) with not-so-distant 2n (from 60 to 72), we assume that the chromosomal evolution 
in Oecomys is more complex than previously thought. In this sense, the use of HME whole chromosome probes 
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MZUSP35534 (2n= 80/FNa= 140) NW 
MZUSP29528 (2n= 80/FNa= 140) S | O. bicolor 
CN285 (2n= 654/FNa= 84) 
LTO05 (2n= 654/FNa= 84) Oepaiies 
ROM114059 
ROM114316 
UU43 ( 2n= 70/FNa= 76) E 
MZUSP29525 (2n= 70/FNa= 74) W 
MZUSP29531 ( 2n= 54/FNa= 54) | O. matogrossensis 
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Figure 6. Bayesian Inference topology based on mitochondrial Cytochrome b and nuclear beta-fibrinogen 
intron 7 concatenated. The numbers above branches indicate posterior probability values for Bayesian 
Inference analysis (only values > 0.90 are shown) and bootstrap values for Maximum Likelihood analysis (only 
values > 60% are shown). Bold numbers indicate the samples from this study. Representatives of different 
lineages within O. bicolor, O. cleberi, O. roberti, O. paricola, and O. catherinae complexes are indicated by the 
letters C (central clade), E (eastern clade), N (northern clade), NW (northwestern clade), S (southern clade), 
W (western clade), and Wm (westernmost clade), following*®. Asterisk (*) denotes specimens identified as 

O. bicolor by*®. Cross (*) denotes karyotype information not obtained from the specimen included in the 
phylogenetic analysis. Sample data are provided in Supplementary Table 1. 




















































































































1 2 3 4 5 6 7 8 9 10 rol 12 13 14 15 16 17 
1 O. auyantepui 2.15 | 2.06 | 2.08 [2.36 | 2.10 | 231 1.87 [2.25 | 214 |236 |231 [224 [232 |243 |223 |229 | 2.26 
2 O. bicolor /O. cleberi 13.60 | 620 | 158 |1.35 | 1.69 | 2.02 1.75 }2.08 | 159 ]184 ]|162 [190 [168 ]1.96 |143 |162 | 1.67 
3 O. catherinae 11.95 | 9.75 | 439 |1.73. | 1.77 | 1.73 1.69 | 1.73 1.99 }2.01 [213 |197 |220 [219 |2.03 |189 | 196 
4 O. concolor 14.16 | 690 | 858 |- 1.58 | 2.41 188 [230 | 195 ]1.99 |207 |195 |2.09 |1.95 |192 |1.75 | 1.80 
5 O. franciscorum / O. 12.57 | 10.84 | 10.54 | 8.33 5.97 | 2.21 1.89 | 2.15 1.89 |2.06 [2.00 /1.73 |189 |1.75 |1.80 |1.73 | 2.02 
mamorae 
6 O. matogrossensis 13.98 | 12.81 8.83 | 13.43 | 14.80 | - 1.94 | 2.18 2.19 | 2.28 2.24 2.37 2.26 2.28 2.41 2.42 2.47 
7 O. paricola 10.67 |11.21 | 920 |9.45 | 10.95 | 11.32 | 3.73 | 1.91 1.86 |220 [2.00 |2.01 |220 |199 |2.03 |1.95 |2.00 
8 O. rex 13.43 112.75 | 8.46 | 12.44 | 12.81 | 11.44 | 10.45 |- 2.21 |231 |226 |2.28 |248 |228 |222 |2.33 | 2.20 
9 O. roberti 12.75 | 9.56 11.36 |9.78 | 11.24 | 12.94 |10.45 | 12.77 | 3.65 |214 [1.78 [2.01 [155 [214 |177 |1.98 |220 
10 | O. rutilus 14.07 | 10.57 | 10.82 |8.96 | 12.81 |11.44 | 12.94 |12.94 |1211 | - 214 [214 |210 [236 |221 |2.09 |202 
11 | O. superans 1420 | 8.71 11.57 |9.45 | 11.69 |11.94 | 10.57 | 11.94 | 8.29 |10.45 | - 237 |218 |222 [2.00 |2.14 | 215 
12 | O. sydandersoni 12.75 |11.07 ]10.95 |8.96 | 9.45 |13.43 | 10.95 | 12.44 |10.45 [10.95 | 12.44 |- 2.08 [1.89 |204 |190 | 2.13 
13 | O. tapajinus 13.23 | 9.39 |12.69 ]9.70 [10.70 |1269 |1269 ]14.18 | 6.72 |10.20 [11.19 ]10.20 Jos0 |213 [1.77 |210 [231 
14 | O. trinitatis 14.65 |12.00 ]13.31 |8.96 | 9.58 |12.44 | 10.32 | 12.44 |11.77 |12.94 |11.94 |8.46 | 10.70 | - 2.26 |2.03 |2.20 
15 | Oecomys sp.1 12.98 | 7.59 |1169 |8.96 | 10.20 |14.43 ]11.82 |11.94 | 846 |1194 [945 |945 |697 |11.94 |- 2.12 | 2.26 
16 | Oecomys sp.2 13.34 | 9.08 | 9.70 |647 | 9.08 |14.43 |10.07 | 12.44 | 9.78 |995 |1045 |846 [10.20 [9.45 |1045 |- 1.94 
17 | Oecomys sp.3 12.75 | 8.40 |10.32 |697 | 11.44 |14.43 |10.45 |11.44 |12.27 |945 [995 ]1045 |11.69 |11.44 11.44 |846 |- 
Outgroup 17.26 |14.57 ]14.05 | 13.43 |16.46 |14.10 |15.05 | 12.60 |14.93 |15.26 | 15.26 | 15.42 |15.09 [15.59 |13.76 | 15.26 | 16.09 








Table 2. Mean genetic p-distances between sequences of Oecomys based on the mitochondrial gene 
Cytochrome b. Intraspecific distances are in bold, and standard deviation values are above diagonal. Values in 
percentage. 
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in representatives of the main lineages of Oecomys provides a more accurate view of chromosomal evolution, 
and associated with detailed phylogeographic studies, is a key factor in understanding the speciation processes 
of this diverse and speciose group of Sigmodontinae rodents. 

By comparing the OAU karyotype with the other Sigmodontinae species investigated by HME whole chromo- 
some probes (Supplementary Table 2), we identified that OAU exhibited the chromosomal signatures proposed 
for the Oecomys genus, namely the fragmentation of HME 1 into three blocks and the syntenic block HME 
(13,22)/21°°. We also detected an exclusive trait for OAU, the syntenic block HME 26/20/18, which is different 
from those described for OCA (HME (9,10)/14/5, 23/19/11 and 26/11) and OPA (HME 4/19); OAU also shared 
and the fragmentation of HME 3 into two blocks with OPA, previously described as exclusive for this species”®. 
Among the eleven chromosomal signatures proposed for the Sigmodontinae subfamily (HME 7/(9,10), 8, 1/12, 
6/21, 11/(16,17), 5/(16,17), 20/(13,22), 15, 19/14/19, 24, and 26)!75°!.69-73, QAU exhibits only five: HME 8, 15, 
24 and 1/12, plus the chromosomal signature HME 19/14/19, which is present as a derived character due to a 
fission that generates the HME 19/14 (OAU 13) and 19 (OAU 12). 


New sex system in Oecomys auyantepui (XX/XY,Y). In the genus Oecomys, sex chromosomes 
variation in size and morphology are frequently related to the addition/deletion of CH, which is common in 
rodents”; the one-armed X chromosome often exhibits CH at the centromeric region, while the bi-armed form 
presents a heterochromatic block in the short arm and the variation in length of the Y is often due to the size of 
the heterochromatic block in the long arm***”*8°°°?, However, the bi-armed X chromosome of O. auyantepui 
from the present study exhibits an euchromatic short arm, with CH at the centromere (Fig. 1b). This indicates 
that events other than the general addition/deletion of CH were involved in this variation, and this is supported 
by the HME 3 hybridization signal in the short arm of OAU Neo-X (Xp) and Y, (see Supplementary Fig. 1), 
validating the emergence of a new sex system (XX/XY,Y,) in O. auyantepui (2n=642654/FNa =84). Regard- 
ing Sigmodontinae rodents, this is the first record of a multiple sex system in Oryzomyini, since this type of sex 
determination system has previously been reported only in representatives of Akodontini (X,X,X,X,/X,X,Y)”, 
Phyllotini (XY, Y,)* and Reithrodontini (XY, Y, and Neo-XY)“ tribes. 

In order to understand the Neo-X origin in O. auyantepui, we constructed a dendrogram that shows a hypo- 
thetical scenario with the chromosomes involved. We made a comparative analysis among the other Sigmo- 
dontinae taxa studied with the same set of probes (Supplementary Table 2) and considered as outgroup the 16 
karyotypes from the genera Akodon, Blarinomys, Necromys, Oxymycterus, Thaptomys (Akodontini), Cerradomys, 
Hylaeamys, and Neacomys (Oryzomyini), while as ingroup we considered the karyotypes of Oecomys catherinae, 
O. paricola and O. auyantepui (Oryzomyini). Except for Cerradomys, all karyotypes from the outgroup showed a 
HME 3 hybridization signal in one large acrocentric chromosome pair. The X was acrocentric, or bi-armed with 
a heterochromatic block in the short arm. Thus, we propose that: (1) the ancestral forms of the HME 3 and X 
chromosomes were medium acrocentrics; (2) O. catherinae maintained the autosomal ancestral form, while the 
X exhibited an addition of CH in the short arm; (3) the HME 3 divided by fission into two blocks of unequal size 
(one large and one small) before the diversification events that led to the O. paricola and O. auyantepui species; 
(4) in O. paricola, the HME 3 small block remained as an independent chromosome, and the X exhibited an 
addition of CH in the short arm; (5) in O. auyantepui, a Robertsonian translocation between the HME 3 small 
block with the acrocentric X formed the Neo-X chromosome (Fig. 7). 

The proposal that intercalary heterochromatic blocks, telomeric repeats and/or rDNA clusters between the 
ancestral X and the translocated autosome served as a boundary that suppresses the X-inactivation process in 
the autosomal portion!®!77*-°° is in accordance with our results, which show a centromeric heterochromatic 
block (Fig. 1b) and a large block of ITS in the Neo-X chromosome of O. auyantepui (Fig. 3). In the rodent Mus 
minutoides (Muridae) immunofluorescence techniques demonstrated patterns of histone modification in three 
types of sex chromosomes (Y, X and a mutant X) that confirmed that the X-inactivation does not spread into 
the translocated autosomal portion’”. This feature is of prime importance and guarantees the viability of this 
multiple sex system. In fact, several studies have described natural populations of vertebrates (e.g., bats, rodents, 
marsupials and ruminants) that possess this type of rearrangement!®!77%30-92.35-3763; in cases where such repetitive 
sequences are absent, deleterious effects such as poor viability and infertility are observed!>'®. The contribution 
of repetitive sequences in the evolution of the X chromosome was also proposed’. The authors microdissected 
the Terricola savii X chromosome into five specific-region probes, hybridized in 20 species of Arvicolinae rodents 
(Myomorpha), and identified multiple intrachromosomal rearrangements, such as centromere shifts, peri- and 
paracentric inversions, which were related to the amplification and distribution of repetitive sequences among 
the Xs of distinct Arvicolinae species’. 

Exceptions from this proposal are documented in the common ancestor of eutherian mammals, since a sex- 
autosome translocation occurred and it may have triggered the separation between Eutheria (placental mammals) 
and Metatheria (marsupials) (166 MYA)°, with no intercalary heterochromatic block found between the X and 
autosomal ancestral segments. Distinct processes from the intercalary heterochromatic block would be involved 
in the regulation of X-autosome viability and in the suppression of deleterious effects’®. 

It is noteworthy that telomeric repeats occur at the ends of chromosomes where they provide stability, while 
ITS are relics of chromosomal fusions that arose during karyotype evolution™. Thus, although we have identified 
other chromosomes resulting from Robertsonian translocations in O. auyantepui (OAU 21-23, 25; Figs. 2 and 
3), the ITS was found only in the Neo-X. An investigation of the effects of these telomeric sequences on gene 
expression, recombination and rearrangements, was made by introducing 800 bp of the telomere repeat (TTA 
GGG) in the adenosine phosphoribosyltransferase (APRT) gene in Chinese hamster ovary cells. The main result 
was that gene rearrangements were greatly increased. This type of chromosome instability is in accordance with 
the proposal that het-ITS (heterochromatic-ITS) seem to be intrinsically prone to breakage“, and that ITSs are 
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Figure 7. Idiograms of the chromosomes involved in the emergence of a new sex-system (XX/XY,Y,) in 
Oecomys auyantepui (OAU; 2n =64265<4/FNa = 84), as assessed based on HME probes™. Idiograms in black (Y 
and Y,) were evaluated only by classical cytogenetics. Arrows show direction of chromosome change. We used 
O. paricola-cytotype A® and O. catherinae-Para”’ as representatives of their respective species. Numbers below 
idiograms correspond to the identification of the chromosomal pair; numbers beside idiograms correspond to 
the HME probes. The bottom box encompasses an idiogram of HME karyotype previously elaborated*® and 
adapted in the present study. Each HME chromosome is shown with a single color, except the pairs (9,10), 
(13,22) and (16,17), which have one color each. (H) Indicates large block of constitutive heterochromatin. 


hotspots for chromosomal rearrangements™. Therefore, the elimination of this sequence during chromosomal 
evolution could be a mechanism that provides karyotypic stability°*®’ and might explain its absence in the rear- 
ranged autosomes of O. auyantepui, while its presence in the Neo-X makes the latter prone to other chromosomal 
rearrangements, despite providing stability against X-inactivation of autosomal segments!®"”. 

We noticed that the other cytotypes found in O. auyantepui (2n=64, 66 and 72) from the Jatapu and Jari 
Rivers (Fig. 8, localities 5 and 4, respectively) exhibited distinct morphologies of the X chromosomes (medium 
submetacentric, large metacentric, and large submetacentric, respectively). Although they are found within a 
simple sex determination system (XX/XY), the X chromosomes had euchromatic short arms and repetitive 
sequences at the centromere, such as an ITS in the karyotypes with 2n = 64, 66°”. Perhaps these differences in 
size and morphology, plus the presence of ITS, could be a clue to a more complex rearrangement of the X chro- 
mosome during its evolution. 

Investigations carried out in some groups with multiple sex systems show that the chromosomal evolution 
in the Neo-X and/or Neo-Y gives rise to other derived systems. This is described in Stenodermatinae bats®, in 
which chromosome painting revealed that a XY, Y, system originated a Neo-XY system, due to a translocation 
between Y, and Y,; this Neo-XY has then diverged into two more derivate systems: in one branch, a fission in 
the Neo-X created a X,X,X,X,/X,X,Y; while in the other branch, a fusion between an autosome and the Neo-Y 
originated a Neo-X,X,X»X,/X,X,Y. In rodents from the genus Reithrodon (Sigmodontinae, Reithrodontini), the 
Uruguay population exhibits a XY,Y, system, while the Brazil population shows a Neo-XY system”; a hypo- 
thetical intermediate Neo-X/Y,Y, formula was the ancestor for the Uruguayan form, while the Brazilian form 
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Figure 8. Map showing the distribution, study area and sampling points for Oecomys auyantepui specimens 
analysed in the present study and from literature. Available karyotypic information is shown with description 
of 2n (diploid number), FNa (autosomal fundamental number), and if is associated with simple (XX/XY) 

or multiple (XX/XY,Y,) sex system. Each symbol indicates the karyotype obtained from the locality; more 

than one symbol indicates that distinct karyotypes were collected at the same locality. Geographic limits of 

O. auyantepui are based on sample points provided in the present study (localities 1-3), and from literature: 
locality 4°’, locality 5°, localities 6-14°°, and from GenBank (localities 15-35). The localities mentioned are: 
Brazil: Para, Obidos, forest fragment 7 km distant from the town’s center (locality 1); Trombetas State Forest 
(locality 2); Grao-Para Ecological Station (locality 3); Jari River (locality 4); Amazonas, Jatapu River (locality 5); 
80 km North of Manaus (locality 6); Amapa, Serra do Navio (locality 7). French Guiana: Trois Sauts (locality 8); 
Paracou (locality 9). Suriname: Brokopondo, Brownsberg Nature Park (locality 10); Sipaliwini, Avanavero Falls, 
Kabalebo River (locality 11). Guyana: Upper Takutu- Upper Essequibo (locality 12), Cuyuni-Mazaruni, Kartabo, 
Cuyuni River (locality 13). Venezuela: Bolivar, Auyan-tepui, South slope, Rio Caroni (locality 14). Samples 
obtained from GenBank (localities 15-35) are detailed in Supplementary Table 1. 


originated through a fusion between the Y, and Y,, that gave rise to the Neo-XY system“. The evolutionary 
process and specific events responsible for the variation in size and morphology of the X chromosomes in the 
O. auyantepui cytotypes (2n=64, 66 and 72)? will be elucidated only after the employment of HME whole 
chromosome probes. 

Taking into consideration the phylogenetic relationships and karyotypic data within O. auyantepui, the COI 
phylogeny is the only one that includes all karyotyped samples for this species from both the present study 
(2n=649650) and the literature (2n = 64, 66 and 72)°*°*. Despite being considered as a valuable tool for high- 
lighting cases requiring systematic attention among small mammal species”, our COI sequences of O. auyantepui 
formed a polytomy with only a few subclades mostly with no satisfactory support (Fig. 4). In fact, specimens of O. 
auyantepui with the multiple sex system grouped in a separated branch only in the concatenated Cytb + FGB-I7 
topologies, but with low support in both BI and ML analyses (Fig. 6). Although our molecular analyses do not 
exhibit a better resolution in the terminal branches, the karyotypic information indicates that we are dealing 
with at least two distinct lineages. We have three karyotypes with a simple sex system XX/XY (2n= 64, 66 and 
72) that could represent variation within a single lineage, while the other lineage corresponds to specimens with 
a multiple sex system XX/XY,Y, (2n=649265<). 

The fact that four potentially karyotypic variants are in a small distribution area (Localities 1-4, Fig. 8) 
indicates that isolating mechanisms are operating, and the rise of a multiple sex determination system may 
be acting as a post-zygotic barrier between these apparently sympatric populations, since the difference in sex 
determination systems will generate aberrations in meiotic synapsis. Taking into consideration the low level of 
genetic divergence within O. auyantepui (mean p-distance 2.15%; Table 2) and the impossibility of interbreeding 
between these two sex systems, it is most likely that a speciation process is already on course. Similar results were 
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observed in Deltamys rodents (Sigmodontinae, Akodontini) where the two distinct sex determination systems 
(XX/XY and X,X,X,X,/X,X,Y) act as a reproductive barrier that is reflected in the phylogenetic divergence of 
11.13-12.14% between the two divergent lineages*”. Despite the low level of genetic divergence in O. auyantepui 
mentioned above, it is worthy of note that the specimen N228 exhibited 5.20% of mean genetic divergence 
from other conspecifics included in the Cytb topology (Fig. 5). Further studies are necessary to evaluate if this 
individual represents a new species, which may be either cytogenetically or morphologically distinct from its 
closely related O. auyantepui individuals. 

Some authors have proposed that changes in the sex determination system can alter behavior patterns and 
contribute to pre-zygotic isolation mechanisms, as described in threespine stickleback fish Gasterosteus aculea- 
tus’®, in which populations with XX/XY and X,X,X,X,/X,X,Y systems exhibit different courtship behavior. In 
the rodent Mus minutoides (Muridae), it was observed that XY females have more reproductive success than 
XX females and are more aggressive and have a stronger bite than XY males’'. Thus, differences in the X chro- 
mosomes of O. auyantepui could lead to behavioral modifications and act in pre-zygotic isolation mechanisms 
between these two taxonomic entities. 

An evaluation of the genetic (Cytb) structure of Oecomys aff. roberti (=O. tapajinus) populations observed 
isolated and stable populations, but with no influence from the mid-Araguaia River opposite banks”. Two sym- 
patric O. paricola “eastern clade” populations (OPA-A and OPA-B) from the Belém area of endemism exhibited 
distinct karyotypes*®. Both works reached similar conclusions for their respective analyzed species, proposing 
that Oecomys taxa can exhibit isolated populations in the absence of strong geographic barriers. 

As discussed above, the fact that four potentially karyotypic variants of O. auyantepui are in a small distribu- 
tion area (Localities 1-4, Fig. 8) could be explained by this type of populational structure for Oecomys taxa®®”. 
In this scenario, the rise of a rearranged chromosomal form within an isolated population could be stablished 
in a few generations, as rodents exhibit a tendency to be organized in demes”’”* with low interbreeding among 
distinct populations. The reproductive barrier would be more intense between the ancestral (XX/XY) and 
derived (XX/XY,Y,) systems, as the hybrid offspring would be infertile’. We propose that we are dealing with 
a case of chromosomal speciation, triggered by the emergence of a new sex system (XX/XY,Y,) in isolated O. 
auyantepui populations. Consequently, O. auyantepui is a species complex with at least two distinct lineages, 
which disagrees with the literature data that recover this taxon as a monotypic species sensu®*. 

A comparative analysis of the G-banding patterns among the Sigmodontinae Neo-X chromosomes of O. 
auyantepui (present work), Reithrodon (Reithrodontini“), and Salinomys delicatus (Phyllotini**) was performed 
and this revealed that the autosomal portion translocated in the Neo-X exhibits similar G-banding patterns in 
the three taxa. This suggests that the same autosomal segment is shared among these distinct lineages. Literature 
data show that distinct genera within groups that have multiple sex systems may also share the same autosome 
in the sex-autosome translocation, as in primates from genera Alouatta” and Aotus’® and in bats from genera 
Artibeus, and Uroderma’’, Chiroderma, Mesophylla, and Vampyriscus®*. However, it has been shown that, in 
rodents from genera Proechimys (Echimyidae) and Nannomys (Muridae), species within a genus can have dif- 
ferent autosomes translocated to the X chromosome”*”». 

We suggest that the employment of HME whole chromosome probes in Reithrodon, Salinomys delicatus and 
Oecomys auyantepui with 2n = 64, 66 and 72 will elucidate the origin of Neo-X chromosomes in Sigmodontinae 
rodents; also, this will shed light on the evolutionary relationships among the four karyotypes of O. auyantepui, 
and clarify if we are dealing with simple (XX/XY) and multiple (XX/XY,Y,) sex determination systems; or with 
derived lineages from the XY,Y, system. 


Material and methods 

Ethics. The specimens were captured using pitfall traps” and kept stress-free with full access to food and water 
until their necessary euthanasia that was performed in accordance with animal welfare guidelines established by 
Brazilian resolution CFMV n.1000/2012. The necessary euthanasia occurred in accordance with animal welfare 
guidelines established by the Animal Ethics Committee (Comité de Etica Animal) from Universidade Federal 
do Para (Permit 68-2015), which also approved all experimental protocols of this research. The captures were 
authorized by the Brazilian Environment Department under license (IBAMA 02047.000384/2007-34). JCP has a 
permanent field permit (number 13248) from “Instituto Chico Mendes de Conservacao da Biodiversidade”. The 
Cytogenetics Laboratory from UFPa has a special permit number 19/2003 from the Ministry of Environment for 
samples transport and 52/2003 for using the samples for research. All methods are reported in accordance with 
ARRIVE guidelines (https://arriveguidelines.org/). 


Samples. We studied the karyotypes of three adult samples of Oecomys auyantepui from distinct locations 
in Obidos municipality, Para state, Brazil. The wildlife samples were collected according to the following: one 
male specimen (UFPAM 2027) was collected in a forest 7 km distant from the town’s center (01° 51’ 15” S, 55° 
32' 53.4" W), one female sample (MPEG 39,927) was collected at the Trombetas State Forest (00° 57’ 45.97” S, 
55° 31’ 20.28” W), and one male sample (MPEG 40,457) was collected at the Grao-Para Ecological Station (00° 
37' 49.01” N, 55° 43’ 42.60” W). The specimens were deposited at the zoological collections of Museu Paraense 
Emilio Goeldii (MPEG) and the Museu de Zoologia da Universidade Federal do Para (MUFPA). Both institu- 
tions are in Belém, Para state, Brazil. 


Cytogenetics. The metaphase chromosomal preparations were obtained from bone marrow extraction”. 
The slides containing chromosomal preparations underwent C-Banding*°, G-Banding*' and FISH® techniques. 
The FISH experiments were performed with human telomeric probes (All Telomere, ONCOR), and with 24 
whole chromosome painting probes from a female of Hylaeamys megacephalus (HME; 2n=54)*; from the 24 
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HME whole chromosome probes, 21 correspond to one chromosome each (including the X chromosome), while 
three corresponded to two pairs of chromosomes each (HME (9,10), (13,22) and (16,17). 


Image capture and analysis. We used an Olympus BX41 microscope and a CCD 1300QDS digital cam- 
era to obtain digital images from G-banded and C-banded karyotypes, which were analyzed using the GenA- 
SIs software v. 7.2.7.34276. The Nikon H550S microscope with a DS-QilMc digital camera captured the FISH 
images, which were analyzed using the Nis-Elements software. The karyotypes were organized according to the 
literature®**. The final images were edited using Adobe Photoshop CS6 software. 


Molecular analysis. We obtained sequences from Oecomys auyantepui samples UFPAM 2027 (Field num- 
ber LTO05) and MPEG 40,457 (field number CN285). We used partial nucleotide sequences of the mitochon- 
drial genes Cytochrome b (Cytb; 801 base pairs) and Cytochrome C Oxidase Subunit I (COI; 657 base pairs), 
and sequence data from nuclear beta-fibrinogen intron 7 (FGB-I7; 727 base pairs). We followed the protocols 
described in® for the extraction, amplification, and sequencing of Cytb and FGB-I7 genes; we also followed the 
same protocols for the COI gene, with the primers Fish F2 (TCGACTAATCATAAAGATATCGGCAC) and Fish 
R2 (ACTTCAGGGTGACCGAAGAATCAGAA)™, The sequences obtained in the present study are available 
on GenBank under accession numbers OM927735, OM927739, OM927737 (CN285); OM927736, OM927740, 
OM927738 (LTO05) (see Supplementary Table 1). 

We also used the sequences available on Genbank (http://www.ncbi.nlm.nih.gov/genbank) from the three 
genes in order to include all COI sequences available in this database for O. auyantepui (48 sequences) plus 
representative COI sequences of Oecomys species (14 sequences), totaling 64 sequences in our COI analysis; all 
Cytb sequences available for O. auyantepui (nine sequences) plus representative Cytb sequences of each Oeco- 
mys species/clade previously recognized**?! (43 sequences), totaling 52 sequences in our Cytb analysis; and 
all FGB-I7 sequences available for O. auyantepui (two sequences) plus representative FGB-I7 sequences of each 
Oecomys species/clade recognized***') (15 sequences), totaling 19 sequences in our concatenated data analysis 
(Cytb + FGB-I7) (Supplementary Table 1). Euryoryzomys nitidus, Hylaeamys megacephalus and Oligoryzomys 
utiaritensis were used as outgroup for Cytb and concatenated analyses, and one sequence of H. megacephalus 
for COI analysis (Supplementary Table 1). 

The alignment and editing of the Cytb, COI and concatenated (Cytb + FGB-I7) sequences were conducted 
using the BioEdit Sequence Alignment Editor program, version 7.0.5.3°° with ClustalW tool. A search for the 
best nucleotide substitution model was made using jModeltest 2.1.6 software** on CIPRES platform®’, which 
selected TIM2 +1+G for Cytb, GIR+1+ G for COI, and TIM2+1+G and TPM2uf+ G for the concatenated 
Cytb and FGB-I7, respectively. 

The phylogenetic reconstructions were made using Bayesian Inference (BI) and Maximum Likelihood (ML) 
methods. BI run in MrBayes 3.2.7** with four chains. The algorithm was based on 50 million generations, sampled 
every 1,000 generations. ML reconstruction was obtained using Garli-2.0° with 1,000 bootstrap replicates and 
majority consensus tree. The tree was displayed and edited in Figtree v. 1.4.3°°. Branches supports were evaluated 
on Bayesian posterior probability in BI and by bootstrap in ML. Cytb genetic divergence values (p-distances) for 
O. auyantepui clades (obtained on Cytb analysis) were calculated with MEGA7”'. 


Map. The map was made using QGIS v. 3.10.7. The shapefiles containing geographic data (elevation, hydrog- 
raphy, and country limits) were obtained from DIVA-GIS”, in the link https://www.diva-gis.org/gdata. Geo- 
graphic limits of Oecomys auyantepui are based on sample points provided in the present study (localities 1-3), 
by Lira® (locality 4), Gomes Junior et al.” (locality 5), Patton et al.*? (localities 6-14), and from GenBank (locali- 
ties 15-35). More detailed information is available in Supplementary Table 1. 


Data availability 

The datasets generated and/or analysed during the current study are available in the GenBank repository (https:// 
www.ncbi.nlm.nih.gov/genbank/). All accession numbers supporting the results reported in this article are found 
in the main text and in the supplementary files. 
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